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LIST OF SYMBOLS 

- wall area 

- cavity length; panel length 

- velocity potential generalized coordinate; see equation III-26 

- cavity width; panel width 

- see equation 111-8^ 

- ambient (equilibrium) speed of sound 

-• cavity depth 

- n th acoustic natural mode 

- acoustic action integral 

- plate action integral 

- see equation III-14 

- neck length of Helmholtz resonator 

- acoustic generalized mass; see equation 'I1I-5 

- structural mass/area 

- normal to surface 

- n** 1 acoustic modal pressure; see equation II I -7 

- pressure 

- m** 1 wall generalized force; see equation III-12 

- vP 1 wall modal coordinate; see equation III-ll 

- n structural generalized mass; see equation III-18 

- acoustic kinetic energy 

- kinetic energy of plate 

- time 
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u 

- fluid particle displacement vector • 


V A 

- acoustic potential energy. 


V P 

strain energy of plate 


V 

- voliaae 


w 

- work 


w 

- wall amplitude 


w 

n 

th 

- n acoustic, mode wall deflection; see 

equation II I -7 

w 

- wall deflection 


x,y,z 

- cartesian coordinates 


z 

- impedance 


£' 

- displacement of fluid 


$ 

- velocity potential 


Y 

- ratio of specific heats 


X 

- see equation VI -33 


$ 

a 

- m th (structural) wall natural mode 


P 

- density 



fch 

- a structural mode (critical) damping 

ratio 

c c . 

*n 

- n acoustic cavity mode (critical) damping ratio 

a 

- frequency 


<,V 

- see equation VI-33 


Subscripts 


o 

- equilibrium 



- perturbation 
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A - absorbent wall; acoustic 

F - flexible wall 

H - Helmholtz resonator 

a,n,r,s - modal numbers 

.R - rigid wall 

Superscripts 

A - acoustic 

a,b,c - cavity a,b,c» respectively 

E - external 

c - acoustic cavity 

* - complementary 
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I INTRODUCTION 

We consider internal sound fields, i.e. sound or pressure variations 
inside of bounded enclosures. We proceed from the simplest geometry to 
the more complex. Specifically the interaction between the (acoustic) 
sound pressure field and the (elastic) flexible wall of an enclosure 
will be considered. One might think of this as the. field of “acoustoe- 
lasticity”. A good introduction to this subject is given in “Sound, 
Structures and their Interaction” by Junger and Feit 1 . Also the author 
has briefly discussed this subject in his book, ’’Aero elasticity of Plates 
and Shells” 2 . 

Fro® the point of view of applications, such problems frequently 
arise when the vibrating walls of a transportation vehicle induce a 
significant internal sound field. The walls themselves may be excited 
by external fluid flows. Cabin noise in various flight vehicles and 
the internal sound field in an automobile are representative examples. 

The first physical model to be considered is a simple one; but one 
which is famous in the acoustic literature. 
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II HELMHOLTZ RESONATOR 

The reader may also wish to consult Kinsler § Frey, Fundamentals 

3 

of Acoustics, pp. 187-207, on this topic. - . 

Consider a "small ,r volume (any linear dimension small compared to 
the acoustic wavelength) so that we may consider all properties constant 
within the volume. For example, consider a box with one vibrating wall 


1 



Ke wish to relate the wall motion to the cavity pressure change and compute 

the coupled structural-acoustic natural frequency. 

From the isen tropic gas law 

p » k p r (1) 

and conservation of mass 

pV « constant (2) 

Consider small changes denoted by ( ) from an equilibrium con- 
dition denoted by Qq. 

p * P q + p 

P * P 0 + P 


V - V + V 
o 


( 3 ) 
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Substituting (3) into (1) and (2) , and linearizing in ( ), 



( 4 ) 




0 


Combining (3) and (4) 



Now 

V * Ad where A is wall area 
o 

V « // wdxdy 

Thus 



(S) 


( 6 ) 


( 7 ) 


This is essentially the “Helmholtz resonator'* approximation. The 
equation of motion of the wall is thus 

3 2 w YP_ 

m // wdxdy (8) 

We have ignored, the stiffness of the wall (which is valid for large 
acoustical stiffness, e.g., small depths) and also ignored acoustical 
inertia (which is valid for low frequencies, e.g., large wavelengths). 
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One can approximately solve (3) using Galerkin's method. 

.Assume 

w - W sin ~ sin (9) 

a t> 

A * ab 

Substitute above into .(8) , multiply by sin sin and integrate 
over the plate area. The result is 


d 2 W W 64 n 

“ dt^ Y p o d 7*- " 0 


( 10 ) 


Substituting W » W e itlJt , we may compute the natural frequency as 


or 


yp o 64 

* ar ^ 


cn) 


Now .the problem we have just analyzed is not the one that Helmholtz 
analyzed. Helmholtz considered a bottle with a neck as in the follow- 



Denote by S the displacement of the fluid in the neck, of the bottle." 
This fluid is considered incompressible and is treated. as a rigid body. 
The equation of motion is 



The natural frequency is 


Y P A 

? o 

&r * 

V p A 
o o 


( 13 ) 


c 2 A 
o 

V A 
o 


It is this problem which is usually termed the "Helmholtz resonator". 

One can also consider the problem of an applied external pressure, 
E 

p , acting on the neck. The equation of. motion is 


p «§ 

o dt z 


Y P o A r E 
§ + p 


(14) 


For the particular case of simple harmonic motion where 


_E- — £ iut 

p n p e 


_ — iut 

i = 5 © 


one may obtain from (14), 


-p i (U 2 ♦ 

pB 0 


y Po A 


Thus an "impedance" may be defined as 



where iti)£ = velocity amplitude and z^ is the impedance of the Helmholtz 
resonator. 

The internal cavity pressure, p c , can be computed from the above 
for a sinusoidal external pressure as 



or » 


1 

P 0 i« 2 V - 1 



(16) 


A similar expression can be obtained for the membrane example 
first considered. The significance of (16) for our purposes is that 
it shows that for 


IV . to 2 < 2 then 

A~ cT 2 


£ 
, E 
P- 


> 1 


Hence for low frequency excitations the internal cavity pressure 
will exceed the external pressure level. 


It is emphasized that the perturbation pressure is assumed not to 

vary within the volume in this model. This is reasonable so long as the 

C 2 IT 

acoustic wavelength, X * o , is much larger than any characteristic 

di 

dimension of the volume. Also the fluid properties are assumed not to 
vary within the neck. 



- 7 - 


HI COUPLED FLU ID- STRUCTURAL MOTION OF AN ACOUSTIC .CAVITY 
WITH A FLEXIBLE AND/OR ABSORBING WALL; 

We consider again the enclosed volume which was considered in 
connection with the Helmholtz resonator problem. 

A direct and concise derivation of the acoustic pressure in a 
closed cavity due to arbitrary motion of the flexible cavity walls 
will be given. The pressure is expanded in terms of the normal modes 
of the rigid-walled cavity. The result, which is valid for any cavity 
geometry, is given in the form of- a set of linear ordinary differential 
equations for the response of each normal mode. For the sake of 
completeness the equations of motion of the flexible wall are also 
derived in terms of structural normal modes and hence the complete 
coupled fluid-structural equations of motion are given. 

Most studies of the acoustic pressure within a cavity due to 
motion of a flexible wall have considered only simple harmonic motion , 
although implicitly the more general case of arbitrary motion can be handled 
via the Fourier integral theorem. However, a very simple result for 
arbitrary wall motions can be obtained directly from the acoustic 
equations using Green’s Theorem without resorting to transform methods. 

The end result of the present method is. an expansion for the acoustic 
pressure within the cavity due to motion of the walls, in terms of the 
normal modes of the cavity with its walls assumed to be rigid. The 
derivation applies equally to all cavity geometries, and as such is a 
generalization of the normal mode expansion ,(or ’’guided wave” expansion. 
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in the terminology used there) derived less directly and for rect- 
angular cavities only in [7J. 

Acoustical Problem 

" Let the cavity occupy a volume V, and be surrounded by a wall surface 
A, of which the portion Ap is flexible* while the remainder \ is rigid. 

If the fluid within the cavity is at rest prior to motion of the wall, 

it 

the fluid pressure p satisfies the familiar wave equation, and associated 
boundary condition: 

6 


iE« - 

3n 


a 2 w 

3t" 


on 


*F 


( 2 ) 


» 0 on A^ 


In these equations p and c are the equilibrium fluid density and acoustic 
velocity within the cavity, and w is the displacement of the flexible 
portion of the wall . in the normal direction n (positive outward) . 

Equation (1) has normal mode solutions F n e lu , n = 0,1,2,... 

1 ‘ 

with the following properties 



( 3 ) 


-r — * 0 on A '(4) 

on 

St /•> 

*This is the perturbation pressure p, of course. The is dropped 
subsequently for convenience. 



r £ n 


(5) 
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~ /,,F F dv = 0 
• V 'V r n 


£ M 


n 


r = n 


A till 

u„ is the n acoustical natural frequency and F^ its related natural 


n 


A A 

mode. Superscript A is to distinguish and from their . structural 
counterparts to be discussed in a later section. 

Note that equation (3) has the solution = 0, F Q » 1. All 

A 

other frequencies o> n , n » 1,2,... are positive and non- zero, however. 

A brief review of acoustical normal mode theory is given in Appendix A. 

The wave equation (1) can be transformed into a set of ordinary 
differential equations in time by using Green's Theorem in the form: 


8F_ 


- v 2 P)dv - / A (p*r - F n §!><* 


( 6 ) 


By defining • 


p = — L-_ 
n P^V 


/ P F n <*v 


(7) 


» n 5 T / “f." 


n Ap 


Ap 


and making use of the fact that p and satisfy equations Cl) and (3), 
and boundary conditions (2) and (4), the following ordinary differential 
equations for the acoustic modes are obtained from (g) : 
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•• i2 

P ♦ « P 
n n n 



( 8 )* 


A dot (’) denotes differentiation with respect to time. The quanti- 
ties P (t) are the coefficients in an acoustical normal mode expansion 
m 

for the pressure: 


P = p ■ l P n F (9)* 

V7 * M A 
n 

Since the normal modes F r satisfy the homogeneous boundary condition 
(4) on the entire wall surface A, the normal derivative of expression 
(9) does not converge uniformly on the flexible portion Ap of the wall 
surface. Expression (9) is suitable, however, for calculating the 
pressure itself throughout the cavity and everywhere on the wall 
surface, including the flexible portion. 

When one of the walls of the, cavity is highly absorbent it is 
usual to model it through a simple point-impedance model where it is 
assumed that 

P * * A *A on A A ( 10 ) 


where the A subscript is used to refer to the absorbent wall 
characteristics, e.g. 

- absorbent wall displacement 
z A - absorbent wall impedance. 


♦Note that for n = 0, we have the Helmholtz resonator result,. L.e. 


F o sl 

W. » /wdA 
M A - ’ 


p o --Vf 

V 

p o c o 2p o ■ - p oV 


/ wdA 


1 


• * P 
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The boundary condition 


3n 


- p 0 w A 


■on A, 
A 


C2) 


still applies, of course. Using (10}, (2) may be written for the ab- 
sorbing wall as 


3£ 

3n 


8 



on A. 


Ca- 


using the above in (6) along with (9), (8) becomes 


P + a> A2 P +A,p_c Z E P C » - A- W 
n n nAOO rnr <F n 

- r V- 

" here c » 5 IL . W* 

A 2 a 




The effect of an absorbent wall is to couple all of the (rigid 

t 

wall) acoustic modes. Of course, as we will see explicitly in what 
follows, the flexible wall also couples all of the acoustic modes 
as well. 

Structural Modal Expansion 

In many technical applications the flexible portion Ap of the 
cavity wall may be a structural element, such as a plate or shell. 
In such cases the wall deflection w is often expressed as a series 
of the form 


„ IB 
m 


( 11 ) 



- 12 - 


in Which the modal functions ¥ are defined over the- region Ap, their - 
properties being determined by' structural considerations.. Using (11) 
a set of structural modal equations describing the wall motion are 
derived. As will be shown subsequently, in these equations the 
cavity acoustic pressure, p, appears in the form of generalized » 
forces 

<W v" C12) 

r 


First, however, the quantities appearing in (8) and (8 A ) 
can be expressed using (7) -and (11) as 


W « Z L a 
n _ na t 


(13) 


where 

L s i / F f dA 
na Ap i n a 

r 

Hence, (8 A ) becomes the acoustic modal equation 

•• \2t * ?• 

P + to P * A.p n c^ l P C * - Z A_ L a 
n n n A O 0 r nr T ran t 

"V ? JT m V~ 

r 


(14) 


(15) 


These can then be solved in conjunction with the companion structural 
equations to be derived in the next section. 


Structural Considerations 

When the structure may be represented by a linear structural 
model the total fluid- structural interaction may be treated in a 
simple way. let the structure be represented by a linear 
(partial, differential) equation. 
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5 w + m |p- - p c - p E (16) 

S is a linear differential operator, representing structural stiffness. 
For example, for an isotropic, flat plate, S = DV 4 , where V 4 is the 
biharmonic operator. The second term on the left-hand side is the 
structural inertia, m being structural mass per unit area. On the 
right-hand side we have two pressure loadings, the first due to the 
cavity acoustics, the second due to some specified "external” 
agent. 

For simplicity, we assume that the structural inodes, ^ , (cf . (11)) 
are normal structural modes satisfying an eigenvalue equation 


• S * B - a d> B 2 » 0 (17) 

and associated orthogonality condition 

/ a 4 dA = M for m » r (18) 

\ a r a 


* 0 for m f r 

is the structural natural frequency and i|> its associated normal 
mode. 


Substituting (11) into (16), using (17), gives 
“» q m ] *■ * - P E 


d 2 a 

*■ ♦ 


(19) 


m 


Multiplying through by and integrating over Ap gives, using (18), 
d 2 q 

tdt 2 ” + “m 2 M a * <£ + <£ 


where (recall (12)) 


<£ 5 /p c 


( 21 ) 
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= - / P E dA 
From (9) and (21-)', 


( 22 ) 


<£ ■ Vo c O 1 P n ^ < 23 > 

” < 

Hence the structural modal equations, (20)', may be rewritten using 
(23) as 


M mC% * •£ ■ p 0=0 *F 2 P n V * 

n t/ 1 - 

n 


(24) 


Summarizing the key relations are (15) and, (24) with determined 
from its definition, (22). (15) and (24) are coupled acoustical-struc- 
tural ordinary differential equations which may be solved by any stand- 
ard method. Since we^ normally are dealing with systems- under forced 

\ 

n 

motion, Q“, the initial conditions on q and P will usually be 
ti m n 

the trivial ones. 

dP 

«a* 3T 


n 


n dt 


at t « 0 


Of course, any physically meaningful initial conditions, could be 
accommodated. 

It should be noted that in applications, the use of normal 
structural modes may not always be convenient and hence direct 
structural coupling may occur. Of course, this will always be true 
for a. nonlinear structural model [7] since the -concept of normal 
modes loses much of its significance there. 

Finally, although the acoustic pressure is obviously the physical 
variable of the greatest interest, some additional insight into the 
nature of the structural-acoustic interaction can be obtained by 
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considering instead the acoustic velocity potential. The two are 
related through Bernouilli's equation 

p * - Pq$ (25) 

If 4> is expanded in its natural modes [they are obviously the 
same modes as for p from the above) 

<p a E a F (26J 

n n 

one can from (9) , recall 


p = pc* Z P n F n 
r on n n 
u n 


M 


n 


(25) and (26), determine that 
=0 


2 P 

cf n =* - 3 


n 




Using (27) in (15) and (24) , one obtains 


VM A 
n 

[a ♦ 
1 n 

A 2 

M* a J 
n n J 

+ .A. 

P c 2 Z 

a C 

A 

00. 

r nr 


( 9 ) 


(27) 


US)* 


M » K * V ■ - | ^ L nm + W 

The underlined terms are the structural -acoustic coupling terms. 

The coefficients are anti-symmetric (to within a multiplicative constant) 
and this characteristic, combined with the first time derivatives of 
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and a^, means that the equations possess gyroscopic coupling. 

Thus it is known that (if -*■ 0 so that there is no absorbing wall) 
the coupled structural - acoustic natural frequencies are real, i.e. 
the system without absorbing walls is an undamped resonator. More- 

g 

over Meirovitch (in a different physical context) has shown how 
one may use conventional matrix eigenvalue methods to determine 
the natural frequencies of gyroscopically coupled systems. 
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IV A VARIATIONAL FORMULATION 

The same results may be obtained through a- 1 variational; formulation. ■ 
9-11 12-13 

Previously Gladwell and later Craggs have considered variational state.- 
ments in the context of finite element representations for the coupled 
structural (wall) - acoustic (cavity) problem. Craggs* work has been 
carried furthest and we use it as a point of departure for our discussiolP. 
Even though we shall obtain the same results- as before, the variational 
statement in its present form is somewhat awkward. Hence the author* s 
preference is for the previous formulation using Green’s theorem for 
the cavity. Others may feel differently, however, or perhaps be in- 
spired to develop a more elegant and/or simpler variational statement. 
Hence the following brief discussion. Following Craggs, consider an 
action integral for the plate 


Ip 2 / [T p - V p ♦ (W g - » A )] dt 
t l 


(28) 


where 


T p - 1/2 // m (|£) 2 dxdy 


Plate Kinetic Energy of Flexible Wall 


or 


v D ff rr 3 2 w. 2 ,a 2 W' 2 , _ i 2 w 3 2 w n,3 2 w 2 , 

Vp ~ + CgyZ-) + 2v + 2(l-v)( T7 ~) Jdxdy 


■p 2 " l "ax" 

E . E 

6W * - // p 6w dxdy 

6W* - // p° • 6w dxdy 

A . . » 

5W » - ]} pg$ 6 w dxdy 


dxdy- 1 

Strain Energy of Wall 

Virtual work due 
to External pressure 


(29) 


Virtual work due 
to Acoustic Cavity 
pressure 


Bernoulli's relation has been used to obtain the last line, i.e. 


f 


- p 0 * 
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Using Hamilton’s Principle on I p (with p E and p c not varied) one obtains 
m w + D V 4 w * ~ p E + p c (30) 


with appropriate boundary conditions on the plate edges involving w 
and its derivatives. 

The action integral for the acoustic system with a prescribed 
motion at the cavity surface is 



< t a- v a - W V dt 


where u - (vector) particle fluid displacement 


(31) 


=* 1/2 /// p Q u ♦ u dxdydz 

V » 1/2 /// dxdydz 

A P 0 C 0 

// w dp dxdy 


Acoustic Kinetic Engery 

(32) 

Acoustic Potential Energy 


is the complementary work done by the virtual forces (^acting 

through the real displacements) at the structural -acoustic' interface. 

Introducing the velocity potential, 

♦ 

U a V$ 

and recalling Bernoulli's relation 
p c = - P Q $ 

the above may be written as 

T a - 1/2 /// P(J (V**V« dxdydz 
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p 

— V. » 1/2 ) 2 dxdydz (33) 

A c Q 

£ **A “ ~ // w Pq 6 ♦ dxdy 

i 

Using Hamilton's Principle with respect to 1^ (with w not varied) 
one obtains 

V 2 $ - 1/ 2 ♦ * 0 (34) 

c 0 

with the boundary condition 


dZ 


W 


on z » 0 (flexible wall) 


(35) 


As the attentive reader will already have noted, otix use of 

cartesian coordinates is a matter of convenience and the results 

readily carry over to cavities and flexible walls of arbitrary shape. 

The necessity for considering both virtual work, SW^, and 

complementary virtual ' 5W^* makes this variational statement appear 

somewhat cumbersome to the author and it is for this reason that 
« 

the earlier approach using Green's Theorem is preferred. 

It is worth noting that an absorbing wall in the cavity can 

i 

also be included, specifically infiW*^. We add a subscript. A, to 
all quantities to denote absorbent. The complementary virtual 
work due to the absorbing wall is 


6W *AA * “ // W A 

Realizing that in Hamilton's Principle the above will appear as 
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:t 2 

t l ' 

one nay integrate by parts to obtain 
Z 2 

f $ W * M dt = /// w A P 0 <5^ dxdy dt 

h 

Further if we model the absorbing wall through a simple point- 
impedance model then 

P ’-*a"a < 3 « 

or . . 

- “o* ■ 2 * »A on a a 


where z^ - absorbent wall impedance 


Thus / 5W* dt * -/// •— $ dxdy dt 

t, ** Z A 


Using the above and proceeding through the formalism of Hamilton's 
Principle one obtains for the boundary condition on the absorbing sur- 
face (e.g. say z = d) 


it * w. 
9z A 


or 

it 

3z 



d 


Modal Equations 

Modal expansions are now made for w and 


( 37 ) 
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w - £ -q^t) * a (x,y) 

<i> * Z F n (x,y,z) 

It is simplest conceptually and generally most efficient 
computationally to -choose the in vacuuo plate modes for iff and the 
rigid wall modes for.F^. Other possibilities exist, of course, but 
they shall not be pursued here. Using the above choice, one has 

T P ■ V2 I M m q^. 
m 

( 38 ) 

V P * V2 < 

m 

where = ffm ip^ 2 dxdy generalized plate mass 

j[ n vacuuo natural plate frequency 
Here we have used the orthogonality of the tp , i.e; 


// m ip tp dxdy * 0 for ui f m 
m n m r n 


Also to obtain V p an integration by parts has been used along with 
the fact that tp^, satisfy the Buler-Lagrange equation 


DV 1 * ip - co 2 m jp 
r a m r m 


0 ■ 


and appropriate boundary equations for on the plate edges;. 
One also has 


s \ - £ \ s %, 


( 39 ) 


where 


- / P E dxdy 
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and iW A - - Ap p I r a n 6^ L 

n m 

where L nm 2 M F n dxdy 
over Ap 

Ap 


( 40 ) 


Ap 5 total area of (flexible) structural wall. 


Using Hamilton’s Principle one obtains, in the usual way, the 
structural modal equations 

* “m 2< U * C - P 0 Z *n L nm ‘ 0 C41) 

n 

Turning now to the determination of the acoustic modal, equations, 
one may construct 


T - 1/2 v > *2 

V. « 1/2 V p 0 2 M A a 2 
A n n 

“o " 

where , 

M n 5 ^/ F m 2 dxd y d z 
V 


(42) 


and 


V = total cavity volume. 




natural (rigid wall) cavity 
frequency 


and * - p n A c J E a 5a L 

A 0 r m m ran 
m n 



t 


1 


5W* ' dt * - A.p 2 f'E E a <5a C dt 
AA A O ' r n nr 

n t 
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where C . = // F F tlxdy 
nr 11 , nr 

on a a 


and = total area of absorbing surface 


In constructing and V^, we have used' the orthogonality of F^, i.e. 

/// P dxdydz = 0 for u> ^ j* w ^ 
r n n r 

X 

and the fact that F^, oj^ satisfy the Euler-Lagrange equation, 

• a 2 

7 2 F ♦ « F * 0 
n n n 



and appropriate boundary conditions on the cavity walls. 
Using Hamilton’s Principle one obtains 

VM ^ [a w ^ a ] - c 2 A_ 
n *• n n o ^F 

+ 

Equations (41) and (43) are 
(15)* which were obtained previously! 


r \ L 

a 


(43) 


ran 


p A c A A, £ a C « 0 
0 0 A r r nr 


the same as equations (24)* and 
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V ACOUSTIC NATURAL .MODES IN MULTIPLY CONNECTED CAVITIES 

For definiteness we consider two rectangular cavities although 
it will be clear that the basic method is applicable to fairly general 
geometries. See sketch. 



X 


Our interest here is in the acoustic natural modes of the two 
cavities connected by an opening of area Ap. In this opening 
there is, in general, a flexible, structural member which deforms 
(as indicated by the dashed line) and thus permits sound transmis- 
sion from one cavity to the next. All other portions of the ca- 
vity walls are rigid. As before the structural member- has a 
displacement-, w, with a modal expansion 

W - s q^t)* (y,z) (1) 

m 

and the acoustical pressures in cavities a and b have modal 
expansions 



( 2 ) 
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The modal equations of notion for the acoustic pressures are then 
(xecal Sections III and IV) 


p a + 
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nm 
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(3) 


and for the structural wal-1 


%]‘ P 0 C 0 A F ‘ P M 


iM 

n (4) 


• - *» i i 

~7 r 

n 

The + sign in the second of (3) and the - sign in the right hand 
side of (4) is a result of the normal to the opening being outward 
from cavity a and inward to cavity b. 

As in Sections III and IV one may introduce a velocity potential 
through Bernoulli’s equation 

(5) 

with a modal expansion 


P a ’ b = - p Q £ a ’ b 


,a,b a,b 

4 ' * I a F 

_ n n 
n 

(3) and (4) then become 


( 6 ) 
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( 8 ) 


g 

(7) and (8) are in standard form for use of Mei-roviteh's algorithm 
to determine the natural frequencies (eigenvalues) and modes (eigen- 
vectors) of the two coupled cavities. 

Of special interest is the case with a pure opening between 
the two cavities,, i.e. one with zero mass and stiffness. In this 
case substantial further simplifications are possible. Thus 

= 0 and (4) becomes (one could deal alternatively with 
(7) and (8)) 


P a L a p b L b 
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To determine the natural frequencies, one assumes 

— iwt 
Sa - Sn e 

p a,b c pfl,b a iait 
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( 10 ) 


Thus from (3) we may determine 
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Substitution of (10) and (11) into (9) gives (before substitution 
one should replace the summation index, m, in (11) by r to avoid 
confusion with the m in (9)) 


l W° 
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r-« 1,2,3,... 


where =11 L a _ L a 


,a n 


nr nm 


M [-u 2 +t»j 

n 1 m 


+ — , I L b L b 
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11 u b r 2 j. Ab 2 , 
M [-liT'+W 1 

n L J 


n - ,m 

The natural frequencies are determined by the condition that 
the determinant of coefficients must vanish 




= 0 


(13) 


(13) is a non-standard eigenvalue problem because of the form 
that a) 2 takes in-Q^, see (12). However, it has one overwhelming 
advantage as compared to the standard (modif i ed .for gyroscopic coupling 8 ) 
eigenvalue approach, embodied .in (7) and (8) : the sire of the matrix 

in (13) will be much smaller than that .in (7) and (8). This is because the 
number of two-dimensional opening modes, ^ , may be much smaller than 
the number of three-dimensional cavity modes, F a , F n b , to achieve a 
given level of accuracy. -This advantage will persist even when the- 
opening is a structural member of finite stiffness and mass, but the 
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mathematics are a bit more cumbersome. 

In closing this section we quote, without detailed derivations, 
two generalizations of the above results. 

Ea 

For an external wall of cavity a whose motion, w , is pres- 
cribed (sinusoidal motion is assumed for simplicity) 

Ea r Ea - 

w * Z q a (14) 

m 

The equation of motion is 


2 Sa ^rm 

B 



(15) 


where 


= ^Ea Z I q£" L‘ 
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2 . Aa 2 , 
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Forced motion will be discussed in much greater depth in Section VI. 

For three cavities, a, b, and c, which are interconnected by 
openings (all other walls rigid) the counterpart of equation (12) ' 
becomes • 
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A ww 

where r - is area of common opening to cavities a and b, etc. 

The generalization to any number of interconnected cavities should 
now be clear. 


The relationship of this analysis to other methods also deserves 


brief mention. Morse and Ingard 6 have considered two coupled cavities 
employing a Green’s Function approach. The final result is the same 
as that obtained here, (12). See Appendix B for details. Alternatively, 
in the context of the variational formulation, as far as computing 
acoustical natural frequencies of multiply connected cavities with 
pure openings and otherwise rigid walls is concerned, the q^ may be 
considered a Lagrange multiplier which enforces the constraint that 
the pressures in cavities a and b must be equal at their common opening. 
Hence the present analysis may be thought of as a component mode synthesis 
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where the components are the individual cavities. The results are 
entirely analogous to those previously obtained in structural vibrations 
except the roles of force and deflection are interchanged* 4 .. In 
structural applications, the Lagrange multipliers are forces of con- 
straint enforcing common deflections where two structural components 
are connected; in acoustical applications, the Lagrange multipliers 
are deflections of constraint which enforce common forces (acoustical 
pressures) where two acoustical components (cavities) are connected. 

With the combined, equivalent single cavity modes determined by 
the method described in this section, it is sufficient to consider . 
a single cavity in the subsequent discussion of forced response in 
Section Y? • 



VI FORCED RESPONSE OF A CAVITY WITH A FLEXIBLE AND/OR ABSORBING WALL 


Basic Model 

Here a single cavity is considered; multiply connected cavities can 
be treated by first determining their combined single cavity natural 
modes as in Section V. 

The total pressure loading on the structural wall is now the 
sum of the internal acoustic cavity pressure and an external (pres- 
cribed) pressure loading. Hence, the structural modal equations 
are 


VS ♦ 2 Vm Sn * W ■ 9m 


CD 


where Q m = I! P¥ n dxdy 


and 


c E 
p + p 


P C * Vo 2 1 V» 




p , the prescribed external pressure, is the new element in the 
theoretical model. may be written 


( 2 ) 
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where Q E m = - // p E dxdy 
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Combining (1) and (3), 

+ 2 ^0^ q a + m m 2 q IQ ]*p^c 0 2 Ap L P n + Q m 
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(5) 
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The cavity modal equations are 
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One can write (5) and (6) in an attractive matrix form by again 
introducing the velocity potential, <j>. Bernoulli's equation is 


P - - P 0 ♦ 


Using an expansion for 
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■ l 4 " F » 


(5) and (6) may be written 
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These equations can be solved numerically by standard methods . 
In the remainder of this section, simplifying assumptions are con- 
sidered which allow useful but approximate analytical solutions 
to be obtained. These should suffice for a .rough estimate and- also 
serve as a guide for more accurate and more elaborate numerical 
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solutions. .In all cases sinusoidal excitation is assumed, i.e. 



where is the excitation frequency. 

Because our model is linear any other time history can be considered 
through superposition of sinusoidal excitation. Random excitation 
may be considered using the sinusoidal excitation results and power 
spectra .methods* 5 . 
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Simplified Models 

One may .anticipate the largest response will. occur when- the 
excitation frequency is near a structural or cavity resonant fre- 
quency . 

(1) Exciting Frequency * Structural Wall Resonant Frequency. 

If a v where s denotes a structural resonant (natural) 
s 

mode, then only q g will be important and all other qj^ may be neg- 
legted. Then, from (6) 


P_ * - A_ L q 
n jF_ ns ^s 

V 

-U) 

s n 

and from. (10) and .(1) 


CIO) 


•• * - 

M [q +2C w.q +u z q ] = 
S in S S s n s s n s J 



Cxi} 


From the right hand side of (11) it is seen that if w n A > 
then the n tn cavity mode gives rise to equivalent mass, while 
if u n < , then the. cavity mode contributes an equivalent 

stiffness. 


Note that the n = 0 cavity mode always contributes equivalent 
stiffness . 

In many practical examples the structural wall resonant.' 
frequencies will be unchanged by the cavity per se, e .g. the sum on 
the Tight hand side of (11) can be neglected. The circumstances 
under which this is not true will be considered in Section VI, part 
(2), Making this assumption for now, from (11) one computes as 
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“ " “s 

q s ■ q s e as) 

M s [2 

From (12), (10) and (2), the cavity pressure may be computed 
25 S 1V 

If « lowest non-zero acoustical resonance, then only U ■ 0 
need be considered in (13) . Hence (13) simplifies to 




Representative numbers are 


Pq V = 1 * «o * 1 . * S “ - 1 
i A. u V 


/ * s dA -2 , / .^ s 2 dA =1/2 


(14)* 
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Thus p c ^ 5 pg, i.e. typically the sound pressure in the cavity will be 

A A 

larger than the external sound pressure., if w < w . For m > u , 

* c c E 

one mist compute p from (13) .and typically p <p ‘ 

(2) .Structural Resonant Frequency Changed by Coupling with Cavity Modes. 

In . most cases where there is a significant change in a structural 
resonant frequency due to the cavity, only the n = 0 cavity mode will 
be important. Hence (11) may be simplified to 

* 2 « s *,v i s j 

" • p 0 c 0 2 V L 0s % * ^ ei “ £t CIS) 

V 

The effect of the cavity (within our approximations) is to modify 
the structural stiffness term, i.e. the total term is now 

M s“s J * “oV *F L 0s 
V. 

Hence (15) may be written in the more compact form 

M s [q s * 2t,»* q s ♦ - q/ (16) 

where the coupled structural -cavity natural frequency, is 

•4 = “s 2 * P 0 C 0 2 *f l Q s - Cl7)* 

MpV . 

. 4 

(17) generalizes an earlier result by Dowell and Voss for structural- 

cavity coupling. Using (17) the results in Section VI, part 1, e.g. (13), 

can now be generalized to include structural-cavity coupling. By 
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examining the ratio of the two terms in (17), one can assess when 
cavity effects on structural resonances are important. 


V‘o 2 a f \ 

m V'oj 2 I 
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If ^ T 2dA 
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Using .our previous estimates for the integrals, the ratio becomes 


p 0 c o 2A f \ l 
a V « s 2 J 


Clearly the lowest frequency structural mode will be most affected 
by the cavity. In extreme cases 5 more than one structural mode may be 
significantly affected. 

Another less frequent case where the structural frequency may be 
changed due to acoustic action is when a. structural mode frequency and 
acoustic mode frequency are close together. Clearly this is an undesire 
able situation, since large sound levels would be anticipated at this 
double resonance. This special case is considered in (4) below. 

(3) Exciting Frequency « Cavity Resonant Frequency. 

E A . 

Whenu -*■ co , then P will be the dominant pressure mode, 
n n 

But (in the absence of absorbent wall damping) if is finite, then 
from (6) one concludes that 
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The structural equation, (5), becomes 
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Solving (19) for gives 
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(20) into (18),. gives 
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Solving (21) for P , 
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( 22 ) 


A 

For the special case where ^ for some m, then the sums 

can be well approximated by a single term and the above relation- 
al 

ship simplifies. Even when <i> n f u m for any m, there may be a pre- 
dominant structural mode in which case the sun could be approximated 
by a single term. First, however, we digress to consider another 
important limiting case . 
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If F n and p are constants over Ap, then a very interesting 
result obtains, namely 


V-P / ■ L nm" F n / V*/ * 

the ratio of sums simplifies and using (2) and (22), 
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Using (4) and (22) in (2) one concludes 
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w = 0 

M 

q.e.d. 


Note that if 
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then the same result holds, namely 
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m 1 J 


Note 


that p c is constant over Ap. Only if there is a single dominant 


structural mode will this simplify further. Now return to (22). 

Single Dominant Structural Mode: 

For a single dominant resonant structural mode, m, £22) becomes 


P - - M A Q E 
n n un 
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c E 


If pH « c F n j on then p” ■ p** on Ap and w = 0 as before. 

Note that, within reason, ^ may be selected as an appropriate linear 
combination of natural modes. Hence a dominant structural mode assumption 
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is quite useful. 

The central conclusion to be drawn from the above is that'-, in 

simple terms, we may think of the cavity , as acting as a vibration 

absorber for the structural wall when the external exciting frequency 

is equal to the cavity resonance frequency. Hence, at most, the internal 
cavity pressure will be equal to the external pressure and indeed, 

“m >:> w n * or m * ^d/or p varies significantly over the flexible 

c E 
wall, then p « p . 


(3a) Effect of Cavity Wall Absorbtion (Damping) . 

In our model the absorbtion characteristics are modeled by an 
impedance, z^. For representative acoustical materials*^ 
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Again we shall assume a single cavity mode dominates. Then an 
equivalent damping ratio of the usual type may be defined 
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Typically 5^ would -be .03 - .3. 

Using an equivalent , (6) becomes 
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As * E , (6a) becomes 
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Using (20) in (26), 
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Solving, 
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Hence effect of cavity damping is to decrease internal pressure field. 
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However, the numerical effect is typically small. 

(4) .Exciting Frequency = Structural Wall and Cavity Resonant 
Frequencies in Close Proximity. 

For simplicity we ignore any absorbing wall'and: consider a single 
dominant structural wall mode and also a single dominant cavity mode 
whose natural frequencies are very near each other. From (9) the 
equations of motion are 
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For a single cavity mode 
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(30), (31) and (32) may be solved to obtain 
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Note that 


* "v PqV 


a Ap 


From (33) the coupled structural-acoustic resonant frequencies 
may be .determined and the associated cavity pressure. For the 


special case of 


ft A _ » 1 
m 


X « 1 (typically X < 1) 


particularly simple results may be obtained. The resonant frequencies 


are then 


» E X 1 ' 2 
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and the cavity pressure at resonance 
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The left hand side of (35) is approximately p /p . Hence for 
typical and X not too small, p C /pE > 1 when a cavity and structural 

wall mode coincide in natural frequencies. 

(5) Convection Effects on External Sound Sources and Internal Sources. 

The results discussed to this point are for a stationary external 
sound source. Here a few brief (cautionary) comments are offered with 
respect to other sound sources. 

If substantial convection of an external sound source is present. 
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hydrodynamic coincidence effects may occur and also radiation 

damping in .the external flow may be important* ^ 

For sound sources inside the cavity itself, the wall absorbtion' 

will play a much more important role than for the external sound . 

X6 

sources studied here . 

(6) Summary of Key Relations from Simplified Models to Determine Cavity 
Pressure. 

Exciting frequency « structural wall resonant frequency. 

P 

For a* + u> ,(13)* or (14)* modified by (17)* as necessary, 
s 

- Exciting frequency = cavity resonant frequency. 

E A 

For uj u * (single dominant structural mode) (25)* 

' n 

Exciting frequency = structural wall and cavity resonant fre- 
quencies in close proximity. 

■ For w E ■+ u = A , (35)*-. 

d n 

When the external excitation frequency is well separated from all 
structural or cavity resonant frequencies, there may not be a" single 
dominant structural and/or cavity mode. However,, some simplification 
may still be possible by neglecting the interaction between structural 
wall and cavity, i.e. one may first determine the structural wall 
motion and then use that result to determine the internal cavity acoustic 
pressure. In this approximation the effect of the cavity sound pressure 
on the structural wall motion is omitted. 

It is also worth emphasizing that for off-resonant exciting fre-r 
quencies, .although the numerics may be more elaborate due to the 
necessity of accounting for multiple modes, the basic theoretical model 
may be more accurate because the uncertainty with respect to structural 
and/or acoustic damping values will be less important. 
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VII NUMERICAL RESULTS AND COMPARISONS WITH EXPERIMENTS 
For a siftgle cavity with a flexible wall and an external sound 
source, the theoretical model has been verified experimentally by 
several authors 4, 5,18,19,20, 21 . 

Here we assess the capability of the model to describe accurately 

the acoustic natural inodes in multiply connected cavities. Once 4 

the combined natural modes of the multiply connected cavities are 

determined and verified experimentally, they may be treated as one 

4- 5 XS 19 21 

single cavity. Hence the earlier work for a single cavity * * * ’ * 

then may be taken as experimental verification for the forced exci- 
tation of multiply connected cavities as well. 

Acoustic Natural Modes in Multiply Connected Cavities . 

?2 

The experimental studies were conducted by Smith" who has con- 
sidered several geometrical arrangements for two acoustic cavities' 
with rigid walls and a partial opening between them, see Fig.l. In Fig. 2 
the lowest found longitudinal resonant frequencies from theory 

(Section V) and experiment are shown vs partition size, y , normalized 

c 

by cavity height, d. As may be seen the agreement is good with the ex- 
ception of the 10* mode for D ^ .4 where a mechanical resonance of 
the loudspeaker used to excite the acoustic modes contaminates the 
data. In practical terms this is unavoidable as both theory and 
experiment show that Q 1Q = .S •+ 0 as the partition is closed, y^ = 

0 ■* .5. In these experiments, c Q » 1117 ft / sec > a ~ d 3 10" a ^d 

the width dimension was 4" to provide two-dimensional conditions in the 

frequency range of interest. The thickness of the partition (assumed 


*10 means there is 1 longitudinal nodeline and 0 lateral notelines, etc. 
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zero in the theoretical calculations) is .5" as is the thickness of 

all external wails. The cavity is constructed from plexiglass. 

In Figure 3 for a constant aperture (opening) size, y . = .5, 

' ' a/ d- 

the effect of aperture location is- studied-. Both’ theory and experiment 
show a small effect. 

A second configuration studied by Smith consists of two’ acoustic 
cavities, one twice the dimensions of the other., with rigid walls and 
a partial opening between them. In Figure 4 the lowest four longitu- 
dinal resonant frequencies from theory and experiment are shown vs. 

normalized partition size, y . . Again the agreement is good. In ‘ 

c/ d 

Figure 5, the longitudinal pressure distributions (along with their 
resonant frequencies) are shown for the first seven (symmetric with re- 
spect to height) acoustical moues with a full opening between cavities. 
The agreement between theory and experiment is very good. 

Altogether one concludes the theoretical results are a faith- 
ful description of the physical model. For more detail the reader 
should consult Reference 22. 

Forced Response of a Cavity with a Flexible Wall 


Experimental Arrangement: 

^For this discussion, Gorman's work^ > ^ is used; however, also 
see Reference 4,5,20 and especially 21. The experimental arrangement 
is shown in Figure 6. The flexible wall panel was a 10” x 20". x 0.05" 
aluminum alloy plate that was bonded onto a rectangular frame consisting 
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of aluminum channel members welded together at their ends. By bonding 
the-.plate to the cavity in this way, a clamped edge boundary condition 
was approximated. A sealed cavity, also 10" x 20'*, was constructed 
beneath the' panel in such a way that the cavity depth could be varied 
in 2 inch increments from 12" to 2" deep. The cavity itself was 
made of 0.5" thick plexiglass and was supported by four plexiglass 
"feet". 

The panel was excited acoustically by a Wolverine LSI 5, 20 watt 
loudspeaker driven by a B § K Beat Frequency Oscillator, type' 1022. 

The external sound field was set at 100 dB at the mid-point on the 
plate surface fox all measurements. By using a single speaker, an 
external field distribution that was variable in space was obtained. 

See Figure 7. 

Initially, there were two basic measurements that were felt 
to be important: the measurement of panel amplitude and the -measurement 

of cavity pressures due to a sinusoidal driving force. As work pro- 
gressed, however, the need arose for two more measurements: the 

measurement of panel and cavity damping ratios. -The reasons for 
thes'e latter measurements will become apparent later. 

Panel Amplitude Measurement: 

The panel motion was measured by the use of a Bently Nevada 
motion pickup. Model 302, that was mounted on an aluminum frame ' 
located above the panel. The frame allowed movement of the pickup 
to any point on the surface of the panel, and also allowed variation 
of the distance between the pickup and the panel. As the panel oscil- 
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lated, the voltage generated by the motion pickup was fed through an 
amplifier and recorded. on a amplitude. ys. frequency plot. .Such a-plot- 
is.sliown in Figure 8, for a cavity depth of 12" . 

t 

For the measurement depicted in Figure 8/ the motion pickup was 
positioned at .the center of the panel . By positioning the pickup in 
this way, one may obtain deflection measurements for the symmetric 
panel modes, i.e,, the modes which have a peak at the panel center, 
but not for the antisymmetric modes, i.e,, those modes with a node 
at the panel center. Hie dominant features' of the plot 

are the three resonant peaks, corresponding to the first, third and 
fifth panel modes, occuring at 113 cps, 210 cps and 410 cps. Modes 
above the fifth mode have an amplitude that is negligible compared to 
the first three symmetric modes. Above 500 cps, the panel 

is essentially motionless. The dominant 

panel response is at the panel fundamental mode, as is expected. 

Cavity Pressure Measurement: 

The sound pressure level within the cavity, when the panel 
has been excited by the loudspeaker, was measured using a B § K 1/4" 
microphone. Type 4136 with a Type 2615 cathode follower with Type UA 
0035 connector. This microphone was installed in holes that were drilled 
in the side of the cavity. As with the panel motion pickup, the voltage 
from the microphone due to excitation of the cavity was recorded on a 
amplitude vs. frequency plot. Such a plot is shown in Figure 9, for 
a cavity depth of 12”. 

\ 

For the measurement depicted in Figure 9, .the microphone was 
positioned , in the hole located at the 3” cavity depth level. Since 
all .theoretical calculations involve the cavity pressure at a point 
just beneath the panel, the 3" depth level was chosen to place the 
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microphone as close to the undersurface of the panel as possible. 

In Figure 9, the cavity pressure is plotted against 
frequency, where this pressure is the difference between the dB 
level inside the cavity and the dB level outside the cavity on the 
upper surface of the panel. Again, the dominant features are the 
three primary resonant peaks occuring at 113 cps, 210 cps and 
SIS cps. The first two resonances correspond to the first and 
third panel modes, and thus indicate that the panel is driving 
the cavity at these frequencies. The resonance at 518 cps is the 
fundamental cavity mode. Note that this mode occurs above the 
frequency at which the panel becomes motionless (500 cps), and thus, 
in effect, the cavity is acting as a rigid cavity with no flexible 
walls. Theoretically, since the panel is motionless, the pressure 
level difference between the external and internal measurements 
should be zero. This is practically the case in Figure 

9, the slight variation from zero due to the very slight motion of 
the panel at this frequency. Note also, that, as with the panel 
amplitude, the greatest cavity response occurs at the panel funda- 
mental frequency, with the second greatest response occurring at the 
cavity fundamental frequency. 

Damping Ratio Measurement: 

The damping ratios of the panel and cavity can be measured 
experimentally in»two different ways: the first, by using the 

"peak method", and the second, by using the "decay method",. The decay 
method was the method used to determine experimentally damping ratios. 
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The actual .measurement o£ damping ratios was as follows: using an. 

oscillograph, a plot of amplitude vs." time. was obtained. As the 
oscillograph was recording, the voltage- to the loudspeaker was cut off, 
thus eliminating the external sound field, and the classic exponential 
decay plot was obtained. 

Figure 10, plots the damping ratio of the panel with a 12" 
deep cavity as a function of frequency. The experimental points 
were measured at the first, third and fifth panel modes for this 
panel - cavity configuration. The curve obtained is consistent 
with a theoretical model, which says that for damping that is propor- 
tional to panel velocity, the damping ratio varies as the inverse of 
the frequency. The slight variation from theory at the higher 
panel frequencies is probably due to coupling between the panel and 
cavity modes. 

Figure 11, plots the damping ratio of the cavity (12" depth} as 
a function of frequency. The experimental points were measured at 
the first, second and third cavity depth modes. The curve obtained 
suggests that the cavity damping ratio at a constant cavity depth 
varies as the inverse of the frequency squared. This result is incon- 
sis tent with Sheshadri 23 , concerning -the damping of Helmholtz resonators, 
and also does not agree with the variation of panel damping with- fre- 
quency at a constant cavity depth. As yet,. no explanation can be given 
for this damping result, other than to conclude that such a result seems 
to .indicate a more complex damping mechanicsm than -that of the panel. 

For the panel, the damping was assumed to be proportional to -the 
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velocity of the panel. For the cavity, such a simple relationship cannot 
be assumed. 

Results : 

Panel Frequency Response 

Since the panel-cavity system being discussed is an integrated 

system, one must investigate the coupling effects between the panel. 

and cavity. It is well known that the cavity effect on the panel will 

become more pronounced as the cavity depth decreases?* * ' In 

Figure 12, the panel deflection, normalized about the mid-point of 

the panel, has been plotted against the panel length, normalized 

with the total panel length, for the fundamental mode of the panel 

backed up by a 12” and a 2" deep cavity. Figure 13, plots a similar 

response, only normalized about the one- third point of the panel, for 

the third panel mode. These graphs depict the effect of intermodal 

coupling at shallow cavity depth?. The panel response at the 12" 

depth for both the first and third panel modes is very close to the 

21 

"in-vacuo" curve given by Dowell and Voss, inter alia, indicating 
-that a 12” deep cavity approximates an infinite cavity fairly effectively 
for the panel size investigated here. 

For a cavity depth„of 2", the mode shapes vary considerably. In 
the case of the panel fundamental mode, the response is more concen- 
trated about the center of the panel than for the 12" depth case. 

The response can be seen to be approaching the shape of the third 
mode, indicated intermodal coupling at this shallow cayity. depth. 

Indeed, if the cavity was made more shallow, the panel, fundamental 
mode would take on the modal pattern of the third mode by having two 
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nodal points^* . For the third panel mode at the 2" cavity depth, 
the change in modal shape is less dramatic; However, there is some 
effect, especially around the mid-point 'of the panel, and this measured 
effect corresponds to. the results of Dowell and Voss. Again, if the 
cavity were made- more shallow, the effect would be more dramatic with • 
the probable elimination of all nodal points, indicating intermodal 
coupling between the first and third panel modes. 

Aside from affecting the mode shapes of the panel, the cavity . 
also affects the natural frequencies of the panel due to the stiffness 
and virtual mass effects discussed previously. In Figure 14, the 
panel fundamental frequency is plotted against cavity depth. Note 
that, as the cavity depth decreases, the panel fundamental frequency 
is an excellent example of the stiffness effect, for, as the cavity 
depth is decreased, the cavity becomes stiffer, thus raising the 
panel natural frequency. Figure 15, plots the panel frequencies vs. 
cavity depth for the third and fifth panel modes. Note that there 
is apparent a slight stiffness effect on the third mode, and the 
lack of any stiffness effect on the fifth mode. In fact, the fifth 
panel mode displays features of the virtual mass effect, for as the 
cavity depth decreases, the natural frequencies- decrease also. This 
fact also indicates that the virtual mass effect is present in sym- 
metric panel modes, 'and not just in the antisymmetric inodes. The 
frequency response for the third and fifth panel modes is consistent ' 
with theory, which states that the stiffness effect becomes less 
for the higher symmetric modes and that the virtual mass effect is pre- 
sent in all modes. 4,5,2i * 
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In Figure 16, a comparison between theory and experiment is made. 

The ratio of panel frequency (modified by coupling with the cavity) to 

"in-vacuo" panel frequency is plotted against panel length to cavity 

depth ratio, a/d. The ."in -vacuo" panel frequencies were computed 

a 

from Warburton's theory, ^ and the panel frequencies* variation with 

21 

cavity depth were computed from Dowell and Voss* theory" which is an 
earlier version of the present analysis. There, is excellent 

agreement between theory and experiment at the large cavity depths, with 
some variation from theory occurring at shallow cavity depths. One 
or two panel modes plus the cavity Helmholtz mode were used in the 
calculation. This agreement between theory and experiment seems to 
indicate that below an a/d of 10, for a panel with similar size, thick- 
ness and material properties as the one investigated here, a one-terra 
panel mode approximation to the panel natural frequencies may be used 
with good accuracy. For a/d > 10, higher term approximations must be 
employed. More generally, for X £ a/d > 10,000 higher term approximations 

must be used to compute panel natural frequencies for a panel length 

S 20 

to width ratio, a/, , of two. ** 

t ’ y 

Cavity Frequency Response 

The cavity also responds in certain characteristic modes. Theo- 
retically, at low frequencies below the cavity fundamental, the cavity 
pressure should be constant over the cavity depth since the cavity is 
responding in its Helmholtz mode. Figure 17, plots the cavity pressure, 
normalized with the external pressure, along the cavity depth at the 
panel fundamental frequency. Apparent is the -fact that the cavity pressure 
distribution is not exactly constant, but is approximately so. 
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Figure 18, plots a similar normalized pressure distribution at 
the cavity fundamental depth frequency. Theory' predicts that the 
cavity should respond in a cosine mode, with the internal and 
external pressures equal at the top and bottom. of the cavity; 
and with the response undergoing a phase shift with cavity depth. 

Note that the experimental pressure distribution in the cavity at 
the cavity depth fundamental follows the theoretical cosine courve 
fairly well except at the deepest part of the cavity. 

Panel' Displacement and Damping Effect 

There are three types of damping that will be referred in this 
discussion: constant damping, frequency damping, and experimental 

damping. Constant damping is the value measured for d = 12" and assumes 

that there is no variation of panel modal damping ratio with cavity 

_ v < 

depth. Frequency damping allows for variation of damping ratio with 

frequency and employs the data measured at a 12" cavity depth for 
various panel resonances. Thus, the only effect changing this type 
of damping is the variation of panel modal frequency with cavity depth. 
Experimental damping is the damping ratio measured for the exact con- 
ditions under investigation. 

/ 

* 

Figure 19, plots the panel damping ratio due to the effect of 
increasing frequency alone, and from experiment. Note that the 
frequency effect forces the panel damping ratios lower, whereas, the 
-experimental damping ratios are higher, as the cavity depth is 
reduced. The reasons for this behavior not readily apparent; however, 
it is felt that, if there were some leakage from the cavity thus 
creating' more losses in the system, the effect would be- to make. 
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the panel respond less ideally, i.e., to increase the damping ratio. 
Experiments show that the cavity pressure increases with -a decrease 
in cavity depth. Thus, if there were a leak present, the effect 
would be greatest at the shallow cavity depths. • This reasoning is 
consistent with the experimental results depicted in Figure 19. 

Kith the panel damping ratios determined, the theory may be 
used to compute a one-term solution for the panel amplitude, for 
three cases: constant damping, frequency damping and experimental 

damping. These results have been plotted in Figure 20, along with 
the measured panel amplitudes. The constant damping 

and frequency damping cases predict the correct panel amplitudes 
only at large cavity depths; however, the experimental damping 
ratios predict the panel amplitudes well throughout the range of 
cavity depths tested. From this, one may conclude, that, given the 
proper damping ratios , a one -term approximation for the panel 

amplitudes is accurate within the range of cavity depths tested. 

18 

A similar result was seen for panel amplitudes for the third 
and fifth panel modes. Again, the experimental damping ratio pre- 
dicted the panel amplitudes more accurately than the frequency 
damping effect. 

Cavity Pressure and Damping Effects 

As with the panel amplitude, the theory may be used to compute 
cavity pressures. Figure 21, plots the- variation of cavity pressure 
with cavity depth for constant damping effect, frequency damping, and 
experimental damping. The exciting external frequency is equal to 
the panel resonant frequency. Recall that the damping ratios used 
in these calculations are the effective panel damping ratios and not 
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those of the cavity; the latter were neglected. Even though 
cavity damping has not been considered,' there is excellent agree- 
ment between experiment- and the experimental damping case. This - 
indicates that, for the fundamental panel mode, and the range of 
cavity depths tested, the use of one-term panel theory, the Helm- 
holtz cavity mode and only panel damping ratios will predict cavity 
pressures accurately. 

Similar results have been obtained for random external pressure 
19 S 

excitation . Pretlove has made measurements of panel natural 

,21 

frequency variation with cavity depth; Guy and Bhattacharya have 
measured cavity pressures and panel natural frequencies. Generally 
good agreement with theory has also been shown in Ref. S, 19 and 21. 
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VIII CONCLUDING REMARKS 

.A comprehensive theoretical model is developed for interior 
sound fields which are created by flexible wall motion resulting 
from exterior sound fields. Including in the model are the. mass, 
stiffness and damping characteristics of the flexible wall as well 
as the mass, stiffness, and damping (due to absorbing interior 
walk) of the acoustic cavity. Full coupling between the wall and 
cavity is permitted although detailed analysis, numerical results 
and experiment suggest that it is the exceptional (though perhaps 
occasionally important} case when the structural w’all dynamics 
characteristics are significantly modified by the cavity. 

Based upon the general theory, an efficient computational 

( 

method is proposed and used to determine acoustic natural fre- 
quencies of multiply connected cavities. 

Simplified formulae are developed for interior sound levels 
in terms of in-vacuuo structural wall and (rigid wall} .acoustic 
cavity natural modes. 

v 

Comparisons of theory with experiment show generally good 
agreement. The principal uncertainty remains the structural and/ 
or cavity damping mechanisms. For external sound excitation, 
cavity damping is demonstrated to be generally unimportant; how- 
ever, it may be of substantial importance for interior sound 
•sources. ‘No systematic experimental data are available for 
cavity damping to assess the validity of the present theoretical 



59 - 


model in this respect. 

The results of the present work, as well as those of Wolf, 
95 

Nefske and Howell using finite element techniques and Howlett 

26 27 

and Morales (based upon earlier work .by Cockburh and Jolly ), 

using modal analysis for a particular wall and cavity geometry, 

suggest that effective analytical models are available subject 

to the uncertainties concerning damping previously discussed. 

For a further entree to the literature one may consult Reference 

25 and additional references cited .therein. as well as the recent 

28 

work by Petyt, Lea and Koopman. 
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APPENDIX A 

ACOUSTIC NORMAL MODES OF A CAVITY WITH. RIGID* WALLS. 

Here we consider a representative, but- simple, calculation of 
acoustical normal modes. 

Normal Modes of a Rectangular Acoustic Cavity 

* 

We shall seek non-trivial (i.e. non-zero) solutions for a 
cavity with rigid walls. As usual the equation of motion is 



with the boundary condition 


|^ = 0 on all (rigid) walls (2) 

where n is coordinate normal to any wall. For a rectangular 
cavity with dimensions and coordinates shown below. 



^2-^0 on x = 0, a 
|£ = 0 on y = 0, b 
8z 


( 2 ) 


0 on z =* 0, d 
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We assume simple harmonic motion 

■p(x,y, z,t) = F(x,y,z) e lut (3) 

and seek the frequencies, w = for which non- trivial solutions 
are possible. These are called the normal mode (or natural) fre- 
quencies. Mathematically speaking they are the eigenvalues of the 
problem. The corresponding solutions for F, F = F m (x,y,z) are the 
normal or natural modes or eigenfunctions. Assuming a separable 

i 

solution for F, 


F(x,y,z) = X(x) Y(y) Z(z) 


(4) 


we may substitute (3) and (4) into (1) to obtain (subscripts denote 
differentiation with respect to that variable) 


X Y 

xx + + 

X Y 



= 0 


Introducing separation constants,' k x . 



(S) becomes 



(5) 


Y 

JQL 

Y 




' ( 6 ) 



= - k 2 (This defines k in terms of 
z z 

a) and vice versa):; 



- 65 - 


Solving (63 

X = A cos k„x + B sin k x 

A A 

r 

Y * C cos k y + D sin k y 

y ■ y 

2 « E cos k z + F sin k z (7) 

z z 

To determine the constants. A, B, etc., k^, etc., we satisfy the 

boundary conditions (2). In terms of X, Y, Z, these are 

X x * 0 on x = 0, a 

Y *0 on y a 0, b 

= 0 on z = 0, d (83 

Using (7) and -the first of these we have 

k x [-A sin + B cos k^xj » 0 x = 0, a 

or k x B » 0 (93 

and k [-A sin k a + B cos k a] » 0 
Thus B = 0 and either 

A a 0 or (103 

sin k x a » 0 -*■ k x a * 0,ir, 2ir, 3ir, ... 

The latter is the alternative which corresponds to non-trivial 
solutions. The corresponding X is 

1 

X = A cos k x = A cos V X , where m = 0,1,2, 3,... (113 

A A 
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Similarly for Y and Z 


kbomir.m =» 0,1,2,... 

y y ' y w w 


Y «* C cos ayny 
b 


( 12 ) 


k d 53 m ir » m = 0,1,2,. 
z z z 


Z * E cos m 7TZ 
z 


Hence from (4) and (6), using (10) , (11), (12), 


F (x,y , 2 ) ~ cos m irx cos m ny cos mvz 
a x y z 

a h . d 


and jo 
c 0 


A -■&)' •(?)'•(?)' 


(13) 


Note F m is only determined to within an arbitrary constant. Also 
we use for shorthand a single modal subscript , m , on F^ and uy. 
Of course* each m corresponds to a single triplet combination of 
integers, ny, m , ny. 

* V 

One of the most interesting and useful properties of the is 
their "orthogonality", i.e. 


/// dxdydz = 0 for m i n 


OH) 


This can be proved most directly using Green's Theorem for and F . 

HI [F V 2 F - F V 2 F J dxdydz 
1 n nit m n J 7 


9F 

// IF. m 


3F n 

n^ 


- F » 3T-J “ 


CIS) 
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3F m 3F 

Now &r - 3?r - 0 


and V 2 f - ' ra 


<i> 2 

- (~> F 


m 'c r m 

(0 ? 

= - • (~) F 
n ^ n 


Hence (15) becomes 


from . (2) 
from (1) 


^n 2 * “m 2) /// F n F m dxd>rdz = 0 
and hence for m f n. 


and /// F^F^dxdydz = 0 

q.e.d. 


( 16 ) 
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APPENDIX B 

STANDING WAVES IN MULTIPLY CONNECTED CAVITIES 

. USING A GREEN’S FUNCTION APPROACH 

An alternative description of the cavity sound field in terms 

of its influence or Green's function is sometimes useful. This will 

be briefly considered here. For further discussion, see Morse and 

6 

Ingard, Chapter 9, pp. 554-564 . 

Green’s Function Determination 

Lef w(x,y,z) = S(x-£) <5(y _T 0 S(z-£) (1) 

Then from II- (7)>W n (t) « 1 F fe.n.O e 1 ** C2) 

II ~i _ I* 

and from II- (8J , = I (iu) 2 O) 

V [(I) 2 - ID 2 ] 

1 U 

Thus from II- (9), p =» l F (5,n,S)F ’ (x,y,z)[-(ia0 2 ]e* ut (4) 

pc/ n 

° ° V “ n *-“ 2 1 

♦ 

The summation 


- O) 2 z F CC,n,?)F (x,y,z) h g 

V n 

M [«o 2 - u> 2 ] 
n n 


may be thought of as a Green's function giving the pressure at x, y, v z 
due to a point harmonic displacement of unit amplitude at 



The Green's function may be used 
to treat two cavities connected 
by a opening. For definiteness 
we consider rectangular cavities, 

although it will be clean the. basic 



method is applicable to fairly general geometries. • 

Consider two volumes or cavities with their respective Green’s, 
functions and a common "wall area", Ap. ' 


G a = - (iu>) 2 Z F a (£,n,5)F n (x,y,z) 

_ a 

V M aTa 2 ^ 

M n t“n - “ 1 


<£> 


I 2 t* n ^ / >- m ^ 


% «.♦ (i«) z 2 F n u (S,n,C) f n (x,y,z) 

“ n. m T 

n ' n 


-rrr 

M u (to - u> 2 ] 


(7) 


There is a + sign on G, because of the sign convention for normal to 
. * o 

Aj.: Also note that only on Ap do both and apply. 

Now 

P a (x,y,z =* 0 on Ap) = // G a (x,y,z = 0 on A F ;£,r>,C * 0) 


p c * 
o o 


over 


A P 


w * 0) d^dn 


(8) 


and 


P b Cx,y,z = 0 on A_)=* // G (x,y,z = 0 on Ap;£,n,£ = 0) 

over Ap 


P.c 


o o 


w a (€,n,5 *0} dcdn 


(9) 


Now further we have on Ap, the physical continuity conditions 


w = w =» w, and p = p, 
a b *a *o 


( 10 ) 


Thus using (8) - (10) 
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This is a homogenuous integral equation for w. It is an eigenvalue 

equation which can be solved in several ways,, e.g. collocation method • 

or modal methods, for the coupled cavities- natural frequencies.- 

■One may also consider further complications due to the wall in 

one cavity moving in a prescribed way. For example at z = - d a 

in room a, say, w * w a d (5,n). Note Ap now includes- both walls, ;=0 

and c * - d . Then 
a 

P a (x,y,z ■ 0 ) a // G a (x,y,z=0;C,n,S=0)w a (5,Ti)dCdn 

Ap 

'-“A 

♦ // C a (K,y,J=0;5,n,C- - d a )w a d (5,n)d?dn (12) 

Ap 

« * - d a 

and from (9), CIO), (12) 

//' (G a -G b ]wd€dri +// G a (x,y,z=0; 5,n,C n “<* a )w a d d£dn = Q (13) 

on A p , on Ap 

?=0 z*0 - 

Given one may solve for w and hence from (8) and (9) , 

determine p in both cavities. In particular if w and are now ex- 
panded in terms of structural inodes, and one uses (6) and (7) in 
(11) and (13), (11) becomes (V-12) and (13) becomes (V-15). These 
same results are derived in a more compact way using the approach 
of Section V.' Hence the present discussion will not be pursued' 
further! 
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